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ABSTRACT 

Models of galaxy formation should be able to predict the properties of clusters of 
galaxies, in particular their gas fractions, metallicities, X-ray luminosity-temperature 
relation, temperature function and mass-deposition-rate function. Fitting these prop- 
erties places important constaints on galaxy formation on all scales. By following gas 
processes in detail, our semi-analytic model (based on that of Nulsen & Fabian 1997) is 
the only such model able to predict all of the above cluster properties. We use realistic 
gas fractions and gas density profiles, and as required by observations we break the 
self-similarity of cluster structure by including supernova heating of intracluster gas, 
the amount of which is indicated by the observed metallicities. We also highlight the 
importance of the mass-deposition-rate function as an independent and very sensitive 
probe of cluster structure. 

Key words: galaxies: clusters: general - galaxies: formation - galaxies: evolution - 
cooling flows - X-rays: galaxies 



1 INTRODUCTION 

Our understanding of galaxy formation through the hierar- 
chical clustering of gas and dark matter (DM) has improved 
greatly over the past decade. The aim of models of galaxy 
formation is to account for the properties of all structure, 
from the first dwarf galaxy to the rich clusters of galax- 
ies that are virialising now. The collapse of the dominant 
dark matter is relatively well understood, so that the ma- 
jor current problem is to model the complex gas processes 
that take place within collapsed and collapsing DM halos. 
Despite their success in accounting for some properties of 
galaxies, many semi-analytical models do not predict cor- 
rectly the properties of groups and clusters of galaxies. In 
particular, these models should predict gas fractions, metal- 
licities, the X-ray luminosity-temperature (Lx — T) relation, 
the temperature function and the mass-deposition-rate (M) 
function of clusters. 

One problem is that the gas fraction in clusters, which 
is observed to be in the range 10-20 per cent (White et 
al. 1993b; White & Fabian 1995; David, Jones & Forman 
1995; White, Jones & Forman 1997; Evrard 1997), is under- 
estimated by those models that have low baryon fractions 
(e.g., most of the values considered by Kauffrnann, White 
& Guiderdoni (1993), Cole et al. (1994), Heyl et al. (1995), 
Kauffrnann & Chariot (1998) and Baugh et al. (1998)). 

In this paper we show that additional heat input to 



the intracluster gas, which in our fiducial model we take to 
be due to supernovae, can consistently account for the ma- 
jor X-ray properties of groups and clusters of galaxies. This 
has a similar effect to 'preheating' as invoked by other au- 
thors (Kaiser 1991; Evrard & Henry 1991; Metzler & Evrard 
1994; Navarro, Frenk & White 1995; Bower 1997; Jones et al. 
1998). A well-known problem addressed by preheating is the 
slope of the Lx — T relation. Models and simulations which 
do not involve preheating (and in general those which inves- 
tigate galaxy properties do not) are generally 'self-similar' 
and give Lx oc T 2 (Kaiser 1986), whereas clusters are ob- 
served to follow L x oc T 3 (David et al. 1993). Models which 
preheat the gas in clusters generally break the self-similarity. 
For example, a constant excess specific energy can have a 
large effect in the shallow potential wells of groups, while 
having little effect on the largest clusters. The excess energy 
raises the temperature as well as the potential energy of the 
gas by pushing it outwards, thereby flattening the gas den- 
sity profile. This lowers Lx and steepens the relationship 
between Lx and T. 

During galaxy formation gas is heated and possibly 
ejected from galaxies by supernovae. It retains most of this 
'excess energy' as thermal and gravitational potential energy 
(and possibly kinetic energy) as it passes through the col- 
lapse hierarchy, eventually becoming part of the intergalactic 
medium in groups and clusters. Gas that is ejected from a 
galaxy is expected to recollapse with the next 'major merger' 
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(which in our model is the doubling, at least, of the mass 
of the halo that the galaxy resides in). This is reasonable 
since ejected gas would encounter the inter-galactic medium 
which subsequently forms part of a larger halo. In our fidu- 
cial model the only source of heating and metal enrichment is 
Type II supernovae (SNe II) associated with star formation 
(Nulsen & Fabian 1995, 1997; hereafter NF95 and NF97). 
There is some evidence that most of the metals in intraclus- 
ter gas came from SNe II rather than SNe la (Mushotzky 
et al. 1996; Nagataki & Sato 1998), but a further reason 
for our assumption is because while both types inject about 
the same amount of energy, SNe la inject about 10 times 
more iron than SNe II. Hence for the iron abundance mea- 
sured in intracluster gas, the associated energy injection is 
maximised by using only SNe II. We assume an energy in- 
jection of 4 x 10 50 erg per supernova (Spitzer 1978) which 
is not radiated locally. The mean mass of iron produced per 
supernova is 0.07 Mq, which leads to 3.7 keV/particle ex- 
cess specific energy in gas with solar iron abundance. By 
choosing the initial metallicity in groups when they are first 
formed to be a suitable value, the iron abundance of clusters 
in our model is Z w 0.3Z© which corresponds to about 1.1 
keV/particle of excess energy. 

The semi-analytical model we use is based on that de- 
scribed in NF97, with an improved cosmological model and 
a more accurate form for dark matter halos. The model pays 
particular attention to gas processes, explicitly including the 
effects of cooling flows in its treatment of hot gas. When hot 
gas has had sufficient time to radiate its thermal energy it 
is assumed to preferentially form low mass stars, effectively 
baryonic dark matter (BDM), as observed in cooling flows in 
clusters. We assume that dark matter halos follow the NFW 
profile (Navarro, Frenk & White 1997; hereafter NFW97). 
Gas in the resulting potential wells is assumed to be isother- 
mal and in hydrostatic equilibrium. The model clusters turn 
out to have gas density profiles closely resembling those in 
observed clusters. 

In this paper we use an open cosmology with fio = 0.3, 
Ho = 50 km s _1 Mpc _1 and a cosmological constant of zero. 
Density fluctuations are obtained from a CDM model with 
primordial spectral index of n = 1, normalised to make erg = 
0.7. This normalisation, which is slightly lower than some 
estimates based on cluster numbers, is discussed in section 3. 
We note that our results would be largely unchanged in a 
flat cosmology with the same fio, since the mass function of 
clusters today would be the same. Likewise, our conclusions 
do not depend on when the excess energies were injected or 
the exact history of the clusters and their progenitors. We 
do however assume that the supernovae injected their energy 
into virialised gas, as opposed to the uncollapsed IGM. 



2 DESCRIPTION OF THE MODEL 

The two main changes from the model described in NF97 are 
the use of an open cosmology and the improvement of the 
DM and gas profiles. The mass of a halo is defined as that 
inside the virial radius, r2oo, and the mean density inside r2oo 
is required to be 200 times the background density of an fi = 
1 universe of the same age. (This is marginally different from 
NFW97, but follows the spherical collapse model strictly.) 
We use the NFW profile (NFW97) which takes the form 



p(r) 



x(l + x) 2 ' 



(1) 



where x = r/r a , and the characteristic density and scale, p s 
and r a , are parameters. For a given mass virialising at some 
time in a given cosmology (which yields r2oo), p(r) has only 
one degree of freedom, which can be fixed by applying the 
algorithm given in the Appendix of NFW97 to obtain the 
concentration, c = 7"20o/?V ps is then fixed since we know the 
mean density inside r2oo- The algorithm used was shown by 
NFW97 to provide reasonable fits of c to simulation results, 
where small halos were found to be more 'concentrated' than 
large ones. In our cosmology the value of c typically ranges 
from 4.5-6.5 for clusters. 

We consider the above profile to describe the total den- 
sity of DM and gas, thus the gravitational potential (j>(r) is 
calculated from it, assuming that the halo is truncated at 
7*200: 



(r) 



ln(l+g) | 1 



1 + c 



(2) 



where a = 4nGp s rg. Since the gas fraction is small, the DM 
also follows the NFW profile to a good approximation. 

In our model the gas is isothermal and in hydrostatic 
equilibrium, giving the gas density profile 



p s (r) oc exp (-tjTjT-^ir)) > 



(3) 



where pmu is the mean mass per particle in the gas, T is its 
temperature and k is Boltzmann's constant. Inserting the 
expression for the potential gives 

p s (r)=A(l + x)^" (4) 

where A is a constant and r\ = pmna./{kT). The above pro- 
file, which we shall refer to as the NFW-gas profile, was 
shown by Makino et al. (1997) to closely approximate the 
conventional /3-model (with /3 = tj/15) . It tends to Ae v in 
the centre and to A for large x. Thus A is, in practice, very 
small and the plateau at large x is not noticeable. The pa- 
rameter a has the same units as velocity dispersion squared, 
thus it is apparent that rj is closely related to the widely- 
used parameter /3, the ratio of the DM kinetic energy to the 
thermal energy of the gas when both are isothermal. 

The NFW-gas profile has been fitted to the surface 
brightness profiles of observed clusters and it models most 
clusters very well (Ettori & Fabian 1998). Values of 77 in the 
largest clusters which have cooling flows (suggesting they are 
in a relaxed state, in contrast to non-cooling flows; Buote & 
Tsai 1996) average about 10.5, with a standard deviation of 
about 1.0 (Ettori, private communication). Direct compar- 
ison shows that the gas profile for 77 = 10.5 closely follows 
the NFW profile from x = 0.4 to 3. The profiles differ in- 
side x — 0.4 due to the fiat core of the NFW-gas profile, 
which implies that the gas fraction increases as a function 
of radius, as observed in real clusters (e.g. White, Jones & 
Forman 1997). This last effect becomes more pronounced 
and spreads to larger radii for smaller values of 77, the major 
cause of which will be heating by supernovae. 

Since 77 gives both the gas temperature and the slope of 
the gas density profile, it is the structural parameter through 
which supernova heating exerts its influence on cluster prop- 
erties. The procedure for setting 77 is as follows. We postulate 
that the total specific energy of the gas (thermal plus gravi- 
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tational) is proportional to the specific gravitational energy 
of the DM, and that any excess specific energy from super- 
nova heating is added to the above gas energy, as its name 
implies, i.e., 



X-ray luminosity vs temperature 



3fcT / P s <t>dV 
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(5) 
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where pdm is the DM density, which follows the NFW pro- 
file, M gas and Mdm are the total gas and DM masses, and 
^excess is the excess specific energy. The volumes of integra- 
tion are spheres of radius r2oo- Given the constant K, it is 
easy to show that i) is only a function of c and E excess . The 
constant of proportionality K is calibrated so that r\ is close 
to the observed value of f0.5 for the largest clusters, when 
^excess = 0. We use the value K = f .2. (Note that r\ varies 
by around 0.5 depending on collapse redshift, through its 
dependence on c.) The reason for matching to the largest 
clusters is that if supernova heating does occur, we expect 
it to have least effect on the largest clusters. As c varies 
little over cluster masses, this method ensures that clusters 
are self-similar when there is no excess energy (i.e., c and rj 
are roughly constant). The trends for c and r\ are to increase 
slightly with decreasing mass, over all masses in our model 
(see below). 

The mass deposition rates of clusters are estimated by 
finding the radius r coo i at which the cooling time is equal 
to the time At to the next major merger (see below) or the 
present, and dividing the gas mass enclosed by At to give 
M. For clusters with cooling flows, i.e., those with significant 
M, the hot gas which cools does not form stars but BDM. 

The Cole & Kaiser (1988) block model is used to simu- 
late the merging history of dark matter halos. The smallest 
masses in the hierarchy are 1.5 x 10 10 Mq and there are 20 
levels of hierarchy so that the total mass in one collapse tree 
is 7.9 x 1O 15 M0. When a 'block' in the hierarchy virialises, 
it merges into one halo both older halos and previously un- 
collapsed material. The new average metallicity and excess 
energy of the gas can therefore be calculated. Since blocks 
double in mass up the hierarchy, each collapse of a block 
roughly corresponds to a major merger. As already men- 
tioned, we set the metallicities in groups to be high without 
considering the implications of this for star formation. Halos 
with masses in the range 15-120 xlO 12 M0 are given initial 
iron abundances of Z = 0.5^0 (This implies that the gas in 
some of the smaller halos is unbound.) As a result the av- 
erage metallicities (for a given mass) range from O.3Z0 for 
clusters of mass 250 x 1O 12 M to 0.25Z e for the largest clus- 
ters. The range of scatter in metallicities is » O.IZq. These 
are in good agreement with the recent ASCA measurements 
by Fukazawa et al. (1998) of iron abundances in 40 nearby 
clusters and the measurements of Allen & Fabian (1998). 

The gas fraction of the clusters in our model lie within 
one per cent of 14 per cent. 



3 RESULTS AND DISCUSSION 

To illustrate the effect of heating from supernovae, we dis- 
play in Fig. ^ the Lx-T distributions before and after the 
inclusion of supernova heating. We superimpose on the dis- 
tributions the power law fit obtained by David et al. (1993), 
for bolometric luminosities. The excess energies bring the 




T (K) 



Figure 1. Contour plots of the luminosity-temperature distri- 
bution, without heating by supernovae (dotted lines) and with 
supernova heating (solid lines). The contours are spaced at equal 
logarithmic intervals. The long straight line is the best fit (for 
bolometric luminosities) taken from David et al. (1993). The dis- 
tribution obtained without heating is roughly parallel to the short 
line segment which follows Lx oc T 2 . In this case the smallest 
clusters are overluminous by an order of magnitude. The excess 
energies bring the smaller clusters into agreement with data. 



smaller clusters into agreement with the observed correla- 
tion, by increasing their temperature, flattening their den- 
sity profiles and thereby decreasing their luminosities. For 
clusters of mass 2.5 x 1O 14 M0, i) decreases to around 9. The 
spread in the Lx-T distributions are significant and arise 
naturally in our simulation from the different formation his- 
tories of the clusters. 

For comparison with Fig. 1, if we halve the excess en- 
ergies injected, the average luminosity of a 2 x 10 7 K cluster 
rises to 5 x 10 43 erg s _1 . Doubling the excess energies causes 
most of the clusters with T < 3x 10 7 K to be unbound, which 
clearly disagrees with data, and roughly halves the luminos- 
ity of 3 x 10 7 K clusters. In future work we will investigate 
whether restricting the heating to the core of clusters can de- 
crease the excess energies required (preliminary results sug- 
gest that this could reduce the total energy requirements by 
about a third). 

In Figs. H and ^ we show the X-ray luminosity function 
(XLF) and temperature function respectively. The curves 
plotted for comparison are the best fit Schechter function 
to the ROSAT Brightest Cluster Sample (BCS) bolometric 
luminosity function, obtained by Ebeling et al. (1997), and 
the best power law fit to the temperature function obtained 
by Edge et al. (1990). The errors in the BCS XLF are very 
small, while the temperature function is currently uncertain 
by factors of around 2 or 3. A problem which follows from 
an incorrect slope of the Lx-T relation is that if the temper- 
ature function fits the data reasonably well, then the XLF is 
too steep compared to data (e.g. Kitayama and Suto 1996). 
Correcting the Lx-T relation thus improves the slope of the 
XLF significantly, though it is still slightly steeper than ob- 
served. 

Estimates of as from cluster abundances in general use a 
theoretically derived mass function and assume some mass- 
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Cluster X-ray luminosity function: z — 0.0 



Mass deposition rate cumulative space density 




Figure 2. The X-ray luminosity function (XLF) from the model 
with supernova heating. The curve is the best fit Schechter func- 
tion to the ROSAT Brightest Cluster Sample (BCS) bolometric 
luminosity function (Ebeling et al. 1997). 

Cluster temperature function; z — 0.0 




Figure 3. The temperature function from the model with su- 
pernova heating. The power law fit to the temperature function 
obtained by Edge et al. (1990) is plotted. 



temperature (M — T) relation in order to compare with an 
observed temperature function. The calibration of the M—T 
relation is thus crucial to the results, especially as the tem- 
perature function has a very steep slope. We find that with 
zero excess energy our model clusters obey the simple esti- 
mate (appropriate for isothermal r~ 2 gas and DM density 
profiles) of 



GMwo _ kT 
2r 2 oo 



(6) 

to within a few per cent, where M200 is the total mass within 
r2oo- Supernova heating raises the temperature slightly for 
small clusters. For a cluster of mass 4.9 x 1O 14 M0 collapsing 
at z = the temperature increases from 3.5 x 10 7 K to 3.8 x 
10 7 K in the model described, close to 10 per cent higher. 
The resulting temperatures in our model are a necessary 




Mass deposition rate (M yr -1 ) 



Figure 4. The smoother curve is the M cumulative function 
from our model with supernova heating, and is compared to the 
function taken from Peres et al. (1997) but with the age of clusters 
changed to 6 Gyr, as explained in the text. We have also excluded 
clusters with T < 2 X 10 7 K; when all clusters are allowed to 
contribute, the function increases below 50 Mq yr -1 , rising to 
6 X 10 -5 Mpc -3 at 10 Mq yr -1 , but such clusters would not be 
included in the observed sample. 



consequence of fitting the Lx — T relation. In a recent paper, 
Eke et al. (1998) used a more elaborate form of Equation ^ 
and included scatter in the M-T relation, which brought 
down the predicted erg. With our CDM shape parameter of 
T = Q Q hexp(-n b (l + 1/fio)) = 0.12 their Fig. 7 predicts a 
value of erg fa 0.75. Since preheating was not included, our 
value of as = 0.7 appears consistent with their results. 

Turning to the mass deposition rates, we display the M 
cumulative space density in Fig. We have excluded clus- 
ters with temperatures less than 2 x 10 7 K since such objects 
would be too faint in the 2-10 keV band to be included in 
the sample used by the two papers below. Our result agrees 
very well with the M cumulative function obtained by Edge, 
Stewart & Fabian (1992) and lies slightly lower than that 
shown in Peres et al. (1997). The systematic uncertainties 
in measurements of M in clusters are at present around a 
factor of 2. In particular the age of a cluster since its last 
major merger is unknown and in determining r coo i the age 
of the universe is generally used. Lowering this age in the 
analysis of Peres et al. (1997) from 13 Gyr to 6 Gyr, which is 
the typical age of a 5 keV cluster in our model, lowers r coo i 
and more than halves the values of M in several clusters 
(Peres, private communication). In Fig. ^| we show the cu- 
mulative function from Peres et al. (1997) recalculated with 
this change in the age of clusters. There is very good agree- 
ment for all but small values of M, which is not a serious 
concern since the clusters and their measurements are both 
sensitive to small perturbations in this range. 

The mass deposition rate of a cluster is a highly sensi- 
tive probe of the gas density inside a radius of ~ 0.1 Mpc 
(the typical r coo i observed in cooling flows). On the other 
hand, the luminosity in all but the most massive cooling 
flows comes from inside a radius ~ r s which is larger than 
the cooling flow region. Therefore Lx and M are indepen- 
dent probes of the gas density profile at different radii, and 
due to their sensitivity they are efficient at eliminating mod- 
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els. For example, while we can adjust the overall density 
inside ~ r s by changing the temperature, the ratio of the 
densities inside the two radii described above have to be cor- 
rect in order to match both the M function and the XLF. 
We found that models which used gas density profiles of the 
form r~ 2/3 (see NF97) overpredicted the M function by an 
order of magnitude or more. The agreement with the avail- 
able data suggests that the NFW-gas profile is much closer 
to the right shape. 

As a further example of the ability of M to discrimi- 
nate between models, without supernova heating we find an 
average value of M « lOOMoyr -1 for cooling flow clusters 
with Lx = 6 x 10 43 erg s _1 which is too high when compared 
with Fig. 2 of Peres et al. (1997). However, heating lowers the 
value to M ~ 30MQyr _1 , which then agrees with their data. 
On the other hand, large clusters with Lx = 10 45 erg s _1 
have the correct M with and without heating. Using a hy- 
drodynamic and N-body simulation, Katz & White (1993) 
modelled a Virgo-like cluster with cooling of gas included, 
and obtained an excessively high value of M ~ 400Moyr~ 1 . 
Our results suggest that non-gravitational heating could si- 
multaneously help solve such a problem. The above results 
argue for flatter gas density profiles in small clusters, in- 
dependently of the Lx — T relation. We consider it very 
promising that by correctly fitting the Lx — T relation, the 
predictions pertaining to M naturally agree with the data. 



4 CONCLUSIONS 

We have presented a galaxy formation model in which the 
X-ray properties of clusters agree well with what is observed. 
The model is based on that described by Nulsen & Fabian 
(1997), and uses observed gas fractions and metallicities. Re- 
alistic gas density profiles were used, which sit in potentials 
given by the NFW profile. In particular, we obtained good 
fits to the Lx — T relation, the X-ray luminosity function, 
the temperature function and the M function, which we find 
to be a very sensitive diagnostic of cluster structure in the 
core. 

To resolve a number of problems, including getting the 
correct slope for the Lx — T relation and avoiding exces- 
sive mass deposition rates for low luminosity clusters, self- 
similarity of clusters has to be broken. In our model this is 
achieved by flattening the gas density profiles of the smaller 
clusters, which occurs naturally when we include heating of 
intracluster gas due to the retention of energy from super- 
novae in all previous collapse stages. We find that the energy 
required agrees well with the supernova energy released in 
producing the observed metallicities if most of this energy 
is retained in the intracluster gas. The result is that the 
above problems are simultaneously resolved once supernova 
heating is included. 

We have highlighted the importance of the mass de- 
position rate and the M function as diagnostics of cluster 
structure. The value of M probes the gas density profile at 
smaller radii than Lx and hence is able to give us indepen- 
dent information. Due to the sensitivity of quantities such 
as Lx and M to the structure and density of the hot gas, 
and the steep slopes of the three distribution functions we 
consider, detailed gas processes can be as important as cos- 
mological considerations in determining what we observe. 
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